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Abstract 

The Kamioka-Korea two detector system is a powerful experimental setup for resolving neutrino 
parameter degeneracies and probing CP violation in neutrino oscillation. In this paper, we study 
sensitivities of this same setup to several nonstandard neutrino physics such as quantum decoher- 
ence, tiny violation of Lorentz symmetry, and nonstandard interactions of neutrinos with matter. 
In most cases, the Kamioka-Korea two-detector setup is more sensitive than the one-detector setup, 
except for the Lorentz symmetry violation with CPT violation, and the nonstandard neutrino inter- 
actions with matter. It can achieve significant improvement on the current bounds on nonstandard 
neutrino physics. 
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I. INTRODUCTION 



Variety of the neutrino experiments, the atmospheric solar [2J, reactor [3|, and ac- 
celerator [ij experiments, have been successful in identifying the neutrino mass induced 
neutrino oscillation as a dominant mechanism for neutrino disappearances. After passing 
through the discovery era, the neutrino physics will enter the epoch of precision study, as 
the CKM phenomenology and a detailed study of CP violation have blossomed in the quark 
sector. The MNS (Maki-Nakagawa-Sakata) matrix elements will be measured with higher 
accuracy, including the CP phase (s), and the neutrino properties such as their interactions 
with matter etc., will be studied in a greater accuracy. During the course of precision stud- 
ies, it will become natural to investigate nonstandard physics related with neutrinos, which 
include flavor changing neutral/charged current interactions 0, @, 0], the effects of quantum 
decoherence jgl, 0, [lO] that may possibly arise due to quantum gravity at short distance scale 
and violation of Lorentz and CPT invariance to name a few. 

It is well known that in history of physics experiments based on interference effect played 
very crucial roles in advancing our understanding of the physical laws. The famous two- 
slit experiment by Young, Michelson-Morley experiment which demonstrated that there is 
no ether, Davidson-Germer experiment on electron diffraction, K° — K° oscillation, etc. 
Likewise, neutrino oscillation experiments may probe another important structure of fun- 
damental physics by observing tiny effects due to CPT violation or quantum decoherence, 
which may be rooted in quantum gravity. Along with the neutral meson systems (K° — K° 
and B° — B°), neutrino oscillations could provide competing and/or complementary infor- 
mations on those exotic effects. See Ref. 14j for a recent review on this subject. 

In a previous work [15j, [l6j we have introduced and explored in detail the physics potential 
of the Kamioka-Korea two detector setting which receive an intense neutrino beam from J- 
PARC. We have demonstrated that the setting is powerful enough to resolve all the ei ght- fold 

19 1, if #13 is in reach of the next generation accelerator 



parameter degeneracy [17j, |18 



and the reactor experiments |22|, |23j. The degeneracy includes the parameters 6 13} S and 
octant of 823, and it is doubled by the ambiguity which arises due to the unknown sign of 
Am^. The detector in Korea plays a decisive role to lift the last one. For related works on 



24, 25, 26, 27 



Kamioka-Korea two detector complex, see, for example, 

The Kamioka-Korea identical two detector setting is a unique apparatus for studying 
nonstandard physics (NSP). As will be elaborated in Sec. [Ill the deviation from the expec- 
tation by the standard mass-induced oscillation can be sensitively probed by comparing 
yields at the intermediate (Kamioka) and the far (Korea) detectors. In this paper, we aim 
at exploring physics potential of the Kamioka-Korea setting in a systematic way. By this 
we mean that we examine several NSP effects in a single framework by concentrating on 
v p, ~ v T subsystem in the standard three-flavor mixing scheme, and focus on disappearance 
measurement. Though limited framework, it will allow us to treat the problem in a coherent 
fashion. 

In analyzing nonstandard physics in this paper, we aim at demonstrating the powerfulness 
of the Kamioka-Korea identical two detector setting, compared to other settings. For this 
purpose, we systematically compare the results obtained with the following three settings: 

• Kamioka-Korea setting: Two identical detectors one at Kamioka and the other in 
Korea each 0.27 Mton 



Kamioka-only setting: A single 0.54 Mton detector at Kamioka 
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• Korea-only setting: A single 0.54 Mton detector at somewhere in Korea. 

Among the cases we have examined Kamioka-Korea setting always gives the best sensi- 
tivities, apart from two exceptions of violation of Lorentz invariance in a CPT violating 
manner, and the nonstandard neutrino interactions with matter. Whereas, the next best 
case is sometimes Kamioka-only or Korea-only settings depending upon the problem. 

This paper is organized as follows. In Sec. [Ill we illustrate how we can probe nonstan- 
dard physics with Kamioka-Korea two detector setting, with a quantum decoherence as an 
example of nonstandard neutrino physics. In Sec. IHIt we describe the statistical method 
which is used in our analyses in the following sections. In Sec. |TV] we discuss quantum 
decoherence. In Sec. |V] we discuss possible violation of Lorentz invariance. In Sec. IVIl we 
discuss non-standard neutrino matter interactions, and the results of study is summarized 
in Sec. IVTU 



II. PROBING NON-STANDARD PHYSICS WITH KAMIOKA-KOREA TWO 
DETECTOR SETTING 

In this section, we describe how we proceed in the following sections. For the purpose 
of illustration, we consider quantum decoherence (QD) as nonstandard neutrino physics. 
In this case, the survival probability (and the 17 M survival probability assuming CPT 
invariance in the presence of QD) is given by 0, To| . 



P(v„ 



P(V„ 



1 - - sin 2 26 
2 



1 _ e ~l(E)L 



COS 



Am 2 L 
2E 



where we consider the case j(E) = j/E (see Sec. IV where we consider also the cases 
of j(E) which has other energy dependences). Then one can calculate the number of 
and events observed at two detectors placed at Kamioka (with baseline of 295 km) and 
Korea (with baseline of 1050 km), using the above survival probability and the neutrino 
beam profiles. For simplicity, let us consider the number of observed neutrino events both 
at Kamioka and Korea, for each energy bin (with 50 MeV width) from E v = 0.2 GeV upto 
E v = 1.4 GeV. In Fig. [TJ, we show the event spectra at detectors located at Kamioka 
and Korea for the pure oscillation 7 = (the left column) and the oscillation plus QD with 
two different QD parameters, 7 = 1 x 10 -4 GeV/km (the middle column) and 7 = 2 x 10 -4 
GeV/km (the right column) 1 . One finds that spectra change for non- vanishing 7. Especially 
we point out that the spectra changes are different between detectors at Kamioka and Korea 
due to the different L/E values at the two positions. 

Assuming the actual data at Kamioka and Korea are given (or well described) by the pure 
oscillation with sin 2 2$ = 1 and Am 2 = 2.5 x 10~ 2 eV 2 , we could claim that 7 = 1 x 10~ 4 
GeV/km (shown in the middle column), for example, would be inconsistent with the data. 
One can make this kind of claim in a more proper and quantitative manner using the y 2 
analysis, which is described in details in the following section. 



1 In order to convert this 7 in unit of GeV/km to 7 defined in Eq. (5), one has to multiply 0.197 x 10" 



18 



3 



sin 2 26=1 .0 , Am 2 =2.5xio _3 (eV 2 ) 



4000 



y =0.0 (GeV/km) 



y =0.0001 (GeV/km) 



y =0.0002 (GeV/km) 




1000 
750 
500 
250 






0.5 


1 


; Korea, 4yr v beam 




n 










0.5 


1 




0.5 1 
Reconstructed ^ (GeV) 



FIG. 1: Event spectra of neutrinos at Kamioka (the top panel) and Korea (the bottom panel) for 
7 = (the left column), 1 x 10~ 4 GeV/km (the middle column), and 7 = 2 x 10~ 4 GeV/km (the 
right column). The hatched areas denote the contributions from non-quasi-elastic events. 



III. ANALYSIS METHOD 

In studying nonstandard physics through neutrino oscillations we restrict ourselves into 
the — v T subsystem, rather than dealing with the full three generation problem. The 
reason is partly technical and partly physics motivated. The technical reason for truncation 
is to simplify analysis in a manner not to spoil the most important features of the problem. 
First of all, — > and — > probabilities do not depend much on the yet unknown 
paramters, namely, $13, CP phase and the sign of Am^. Second, in most cases the earth 
matter effect is a sub-leading effect in the — u T subsystem so that by restricting to the 
subsystem we need not to worry about complications due to the matter effect. Moreover, 
many of the foregoing analyses were carried out under the truncated framework. By working 
under the same approximation the comparison between ours and the existing results becomes 
much simpler. Thus, the Am 2 which will appear in neutrino oscillation probabilities in the 
following sections is meant to be Am 2 = Am| 2 = m 2 — m|. 6 will be a mixing angle in the 
unitary matrix diagonalizing the Hamiltonian in the — v T subsystem, which is essentially 
equal to #23 . 
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A. Method of statistical analysis 



In order to understand the sensitivity of the experiment with the two detector system at 
295 km (Kamioka) and 1050 km (Korea), we carry out a \ 2 analysis. In the present analysis, 
we only included and I7 M disappearance channels. In short, the definition of the statistical 
procedure is similar to the one used in [l6j excluding the electron events. The assumption on 
the experimental setting is also identical to that of the best performance setting identified 
in Ref. [IBj]. Namely, 0.27 Mton fiducial masses for the intermediate site (Kamioka, 295 km) 
and the far site (Korea, 1050 km). For the reference, we also consider 0.54 Mton detector 
for Kamioka or Korea only. The neutrino beam is assumed to be 2.5 degree off-axis one 
produced by the upgraded J-PARC 4 MW proton beam. It is assumed that the experiment 
will continue for 8 years with 4 years of neutrino and 4 years of anti-neutrino runs. 



in 



We use various numbers and distributions available from references related to T2K [28 
which many of the numbers are updated after the original proposal [loj . Here, we summarize 
the main assumptions and the methods used in the \ 2 analysis. We use the reconstructed 
neutrino energy for single- Cherenkov-ring muon events. The resolution in the reconstructed 
neutrino energy is 80 MeV for quasi-elastic (QE) events. We take Am 2 = 2.5 x 10~ 3 eV 2 
and sin 2 2#23 = 1.0 for our reference value. However, whenever we expect that there is a 
correlation between the expected sensitivity to new physics and the oscillation parameters, 
we scan Am^ between 2.0 and 3.0xlO _3 eV 2 and sin 2 2^23 between 0.9 and 1.0. The shape 
of the energy spectrum for the anti-neutrino beam is assumed to be identical to that of 
the neutrino beam. The event rate for the anti-neutrino beam in the absence of neutrino 
oscillations is smaller by a factor of 3.4 due mostly to the lower neutrino interaction cross 
sections and partly to the slightly lower flux. In addition, the contamination of the wrong 
sign muon events is higher in the anti-neutrino beam. 

We stress that in the present setting the detectors located in Kamioka and in Korea 
are not only identical but also receive neutrino beams with essentially the same energy 
distribution (due to the same off-axis angle of 2.5 degree) in the absence of oscillations. 
However, it was realized recently that, due to a non-circular shape of the decay pipe of the 
J-PARC neutrino beam line, the flux energy spectra viewed at detectors in Kamioka and 
in Korea are expected to be slightly different even at the same off-axis angle, especially 
in the high-energy tail of the spectrum [23|. The possible difference between fluxes in the 
intermediate and the far detectors is taken into account as a systematic error in the present 
analysis. 

We compute neutrino oscillation probabilities by numerically integrating neutrino evolu- 
tion equation under the constant density approximation. The average density is assumed to 
be 2.3 and 2.8 g/cm 3 for the matter along the beam line between the production target and 
Kamioka and between the target and Korea, respectively [3]. We assume that the number 
of electron with respect to that of nucleons to be 0.5 to convert the matter density to the 
electron number density. 

The statistical significance of the measurement considered in this paper was estimated 
by using the following definition of x 2 : 



20 



x 2 



E E J ™ ' +E If.) . ( 2 ) 



k=l \i=l 1 J j=l v ■> 
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where 



N(fi)T P = K° n ~ QE ■ (1 + E /(A*)} • ^ + N? E -(1+ £ /(A*)} " e,-) • (3) 

.7=1,3,4 3=1,2,4 

In Eq. (ED, iV(/i)° bs is the number single-ring muon events to be observed for the given 
(oscillation) parameter set, and iV(yu)° xp is the expected number of events for the assumed 
parameters in the x 2 analysis, k = 1,2,3 and 4 correspond to the four combinations of the 
detectors in Kamioka and in Korea with the neutrino and anti-neutrino beams, respectively. 
The index % represents the reconstructed neutrino energy bin for muons. The energy range for 
the muon events covers from 200 to 1200 MeV. Each energy bin has 50 MeV width. Oi denotes 
the statistical uncertainties in the expected data. The second term in the x 2 definition 
collects the contributions from variables which parameterize the systematic uncertainties in 
the expected number of signal and background events. 

jy-non-QE arg ^ e nU mber of non-quasi-elastic muon events for the i th bin whereas N® E are 
the number of quasi-elastic muon events. We treat the non-quasi-elastic and quasi-elastic 
muon events separately, since the neutrino energy cannot be properly reconstructed for non- 
quasi-elastic events. Both N® E and j\r non_ Q E depend on neutrino (oscillation) parameters 
but in a different way, namely, the former (latter) being affected in direct (indirect) manner 
and hence dependence is strong (weak) as we can see from the solid histgram 
jynon-QE^ ^g hatched region (j\r non- Q E ) j n pig [jj The kgy to high sensitivity to NSP is 
the different oscillation parameter dependence between Kamioka and the Korean detectors 
due to different baselines. The uncertainties in j\r non_ Q E an( j _/Vj9 E are represented by 4 
parameters €j (j = 1 to 4). 

During the fit, the values of A^(/i)^ xp are recalculated for each choice of the (oscillation) 
parameters which are varied freely to minimize x 2 , and so are the systematic error parameters 
€j. The parameter /(/•«)] represents the fractional change in the predicted event rate in the 
i th bin due to a variation of the parameter 6j. We assume that the experiment is equipped 
with a near detector which measures the un-oscillated muon spectrum. The uncertainties in 
the absolute normalization of events are assumed to be 5% (5"i=0.05). The functional form 
of f(fj,) l 2 = (E u (rec) — 800 MeV)/800 MeV is used to define the uncertainty in the spectrum 
shape for quasi-elastic muon events (er 2 =0.05) fl6| - 

The uncertainty in the separation of quasi-elastic and non-quasi-elastic interactions in 
the muon events is assumed to be 20% (03=0. 20). In addition, for the number of events 
in Korea, the possible flux difference between Kamioka and Korea is taken into account in 



The predicted flux difference 29] is simply assumed to be the 1 a uncertainty in the 
flux difference (04). 

Finally, in this work, the sensitivity at 90% (99 %) confidence level (CL) is defined by 
A-X 2 = Xmin( osc - + nonstandard physics) — Xmi n ( osc - 

) > 2.71 (6.63), 

corresponding to the one degree of freedom. Similarly, the criterion for the two degrees of 
freedom is Ax 2 > 4.61 (9.21). 



IV. QUANTUM DECOHERENCE (QD) 

Study of "quantum decoherence" is based on a hypothesis that somehow there may be 
a loss of coherence due to environmental effect or quantum gravity and space-time foam, 
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etc |8j. We do not discuss the origin of decoherence in this paper, but concentrate on how 
this effect can be probed by the Kamioka-Korea setting (this statement also applies to other 
nonstandard physics considered in this paper). For previous analyses of decoherence in 
neutrino experiments, see e.g., 0, El, 32, 33]. As discussed in Sec. H] we consider the 
V/j, — v T two-flavor system. Since the matter effect is a sub-leading effect in this channel we 
employ vacuum oscillation approximation in this section. The two-level system in vacuum 
in the presence of quantum decoherence can be solved to give the survival probability 



1 - - sin 2 29 
2 



e-^ L cos 



Am 2 L 
2E 



(4) 



where 'j(E) > is the parameter which controls the strength of decoherence effect. Notice 
that the conventional two-flavor oscillation formula is reproduced in the limit ^(E) — > 0. 
Since the total probability is still conserved in the presence of QD, the relation P(z/ M — > v T ) = 



1 - P(u, 



holds. A similar expression holds for the anti-neutrino survival probability 



with possible different decoherence index 7(E). We assume CPT invariance in this section 
so that ■y(E) = j(E); Then, the survival probability is the same as in (jlj) follows. 

Unfortunately, nothing is known about the energy dependence of j(E). Therefore, we 
examine, following [||, several typical cases of energy dependence of j(E): 



l{E) = 7 



E 



GcV 



(with n = 0,2,-1) 



(5) 



In this convention, the overall constant 7 has a dimension of energy or (length) 1 , irrespective 
of the values of the exponent n. We will use 7 in GeV unit in this section. In the following 



three subsections, we analyze three different energy dependences, 



n 



0, —1, 2 one by one. 



A. Case of 4 dependence of 'y(E) 

First, we examine the case with ■= dependence of 'j(E). In Fig. El we show the correlations 
between 7 and sin 2 29 at three experimental setups for the case where the input values are 
7 = 0, Am 2 = 2.5 x 10~ 3 eV 2 and sin 2 29 = 0.96. We immediately find that there are 
strong correlations between sin 2 29 and 7 for the Kamioka-only and Korea-only setups. We 
also note that the slope of the correlation for the Kamioka-only setup is different from that 
for the Korea-only setup. Therefore the Kamioka-Korea setup can give a stronger bound 
than each experimental setup. This advantage can be seen in Fig. [31 where we present the 
sensitivity regions of 7 as a function of sin 2 26 l (left panel) and Am 2 (right panel). Note 
that the sensitivity to 7 in the Kamioka-Korea setting is better than the Korea-only and 
the Kamioka-only settings by a factor greater than 3 and 6, respectively. 



B. Case of energy independent 7(E) 

Next, we examine the case of energy independent 7, "f(E) = 7 = constant. Repeating 
the same procedure as before, we find that the sensitivity to 7 in the Kamioka-Korea setting 
is better than the Korea-only and the Kamioka-only settings by a factor greater than 3 and 
8, respectively, as summarized in Table [0 
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FIG. 2: The correlations between 7 and sin 2 20 for three experimental setups we consider: 
Kamioka-only, Korea-only and Kamiok-Korea. Blue, black and red curves represent the contours 
for 68% , 90% and 99% CL for two degrees of freedom. Input values are 7 = 0, Am 2 = 2.5 x 10" 3 
eV 2 and sin 2 20 = 0.96. 



C. Case of E 2 dependence of j(E) 

Finally, we examine the case with E 2 dependence of 7. The qualitative features of the 
sensitivities are similar to those of the previous two cases. The result is that the sensitivity 
to 7 in the Kamioka-Korea setting is better than the Korea-only and the Kamioka-only 
settings by a factor greater than 3 and 5, respectively. 
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FIG. 3: The sensitivity to 7 as a function of the true (input) value of sin 2 29 = sin 2 2023 (left panel) 
and Am 2 = Am 2 2 (right panel) for the case of 1/E dependence of 'y(E). The red solid lines are for 
Kamioka-Korea setting with each 0.27 Mton detector, while the dashed black (dotted blue) lines 
are for Kamioka (Korea) only setting with 0.54 Mton detector. The thick and the thin lines are 
for 99% and 90% CL, respectively. 4 years of neutrino plus 4 years of anti-neutrino running are 
assumed. In obtaining the results shown in the left and right panel, the input value of Am 2 2 is 
taken as +2.5 x 10~ 3 eV 2 (with positive sign indicating the normal mass hierarchy) and that of 
sin 2 26*23 as 0.96, respectively. 



D. Comparison between sensitivities of Kamioka-Korea setting and the existing 
bound on 7 

In Table HJ we list, for the purpose of comparison, the upper bounds on 7 at 90% CL 
obtained by analyzing the atmospheric neutrino data in j^] 2 . We summarize in the table the 
bounds on 7 at 90% CL achievable by the Kamioka-only, the Korea-only and the Kamioka- 
Korea settings. We use the same ansatz as in j^] for parameterizing the energy dependence 
of 7(E), n = 0,-1 and 2. 

In the case of dependence of ^{E), all three settings can improve the current bound 
almost by two orders of magnitude. Note that the best case (Kamioka-Korea) is a factor of 
6 better than the Kamioka-only case. This case demonstrates clearly that the two-detector 
setup is more powerful than the Kamioka-only setup. We notice that in the case of energy 
independence of ^(E), Kamioka-Korea two detector setting can improve the current bound 



2 We do not quote the bounds on 7 obtained from solar and KamLAND neutrinos, since they are derived 
from v e 



v.. 



34j . In this case, the neutrino energy is quite low, so that the constraint for n = — 1 



becomes quite strong. 
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Ansats for j(E) 


Curent bound (GeV) 


Kamioka-only (GeV) 


Korea-only (GeV) 


Kamioka-Korea (GeV) 


j(E) = 7 (const.) 
y(E) = 7 /£(GeV) 
j(E) = j(E(GeV)) 2 


< 3.5 x 10- 23 

< 2.0 x 1(T 21 

< 0.9 x 10" 27 


< 8.7 x 10~ 23 

< 4.0 x 10" 23 

< 9.2 x 10" 23 


< 3.2 x 10" 23 

< 2.0 x 10" 23 

< 6.0 x 10" 23 


< 1.1 x 10~ 23 

< 0.7 x 10" 23 

< 1.7 x 10" 23 



TABLE I: Presented are the upper bounds on decoherence parameters 7 defined in (|5j) for three 
possible values of n. The current bounds are based on [9] and are at 90% CL. The sensitivities 
obtained by this study are also at 90% CL and correspond to the true values of the parameters 
Am 2 = 2.5 x 10~ 3 eV 2 and sin 2 2<9 23 = 0.96. 



by a factor of ~ 3. 

In the case of E 2 dependence of ^(E) the situation is completely reversed; The bound 
imposed by the atmospheric neutrino data surpasses those of our three settings by almost 
~4 orders of magnitude. Because the spectrum of atmospheric neutrinos spans a wide 
range of energy which extends to 100-1000 GeV, it gives much tighter constraints on the 
decoherence parameter for quadratic energy dependence of ^y(E). In a sense, the current 
Super-Kamiokande experiment is already a powerful neutrino spectroscope with a very wide 
energy range, and could be sensitive to nonstandard neutrino physics that may affect higher 
energy neutrinos such as QD with j(E) ~ E 2 or Lorentz symmetry violation (see Sec. V for 
more details). 



V. VIOLATION OF LORENTZ INVARIANCE 

In the presence of Lorentz symmetry violation by a tiny amount, neutrinos can have 



both velocity mixings and the mass mixings, both are CPT conserving [ll|]. Also there 



could be CPT- violating interactions in general [111, ll2|, 113] . Then, the energy of neutrinos 



with definite momentum in ultra-relativistic regime can be written as 



mm) m 2 



where m 2 , c, and b are 3x3 hermitian matrices, and the three terms represent, in order, the 



effects of velocity mixing, mass mixing, and CPT violation 12]. The energies of neutrinos 
are eigenvalues of and the eigenvectors give the "mass eigenstates" . Notice that while c 
is dimensionless quantity, b has dimension of energy. For brevity, we will use GeV unit for 
b. 

Within the framework just defined above it was shown by Coleman and Glashow 12| 
that the disappearance probability in vacuum can be written as 

P(v„ -> = 1 - sin 2 26 sin 2 (AL/4) (7) 

where the "mixing anlge" G and the phase factor A depend upon seven parameters apart 
from energy E: 

Asin2G = Am 2 sin2^ m / J E + 256e 4 ' ? sin2^ + 25ce ir '' J Esin2^ c , 

Acos2Q = Am 2 cos2^ m / J E + 256cos2^ + 25c J Ecos2^ c . (8) 
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As easily guessed what is relevant in neutrino oscillation is the difference in mass squared, 
and 6 and c between two mass eigenstates, 5b = 62 — 61 and 5c = c% — ci, where Cf=i ( 2 
and 64=1,2 are the eigenvalues of the matrix c and 6. The angles # m , 9b and 6* c appear in 
the unitary matrices which diagonalize the matrices m 2 , b, and c, respectively. There are 
also two phases 77 and rj that cannot be rotated away by field redefinition. We work in the 
convention in which cos 29 m and cos 29b are positive, and Am 2 = m 2 , — mf, 5b = b 2 — bi and 
5c = C2 — Ci can have either signs. 

The survival probability for the anti-neutrino is obtained by the following substitution: 

5c -> <5c, 56 -> -56 (9) 

The difference in the sign changes signify the CPT conserving vs. CPT violating nature of 
c and b terms. 

The two-flavor oscillation given by ([7]) and (jSJ) is a too complicated system for full analysis. 
Therefore, we make some simplifications in our analysis. We restrict ourselves into the case 
9 m = 9b = 9 C = 9 and rj = rj = 0, for which one recovers the case treated in [35j: 



P(v» -»• f M ) = 1 - sin 2 29 sin z 



L 



Am 2 5b 5cE 



AE 2 



(10) 



which still depends on 4 parameters, 9, Am 2 , 5b and 5c. One has a similar expression for V M 
with 5b — > —<56. As pointed out in [36j, the analysis for violation of Lorentz invariance with 5c 
term is equivalent to testing the equivalence principle [37]. The oscillation probability in (flUj) 
looks like the one for conventional neutrino oscillations due to Am 2 , with small corrections 
due to the Lorentz symmetry violating 5b and 5c terms. In this sense, it may be the most 
interesting case to examine as a typical example with the Lorentz symmetry violation. Note 
that the sign of 5b and 5c can have different effects on the survival probabilities, so that the 
bounds on 5b and 5c could depend on their signs, although we will find that the difference 
is rather small. 

For ease of analysis and simplicity of presentation, we further restrict our analysis to the 
case of either 5b = and 5c 7^ 0, or 5b 7^ and 5c = 0. Notice that the former is CPT 
conserving while the latter is CPT violating. 



A. Case with 5b = and 5c / (CPT conserving) 

We first examine violation of Lorentz invariance in the case of 5b = and 5c 7^ 0. In Fig. H] 
we present the region of allowed values of 5c as a function of sin 2 29 (left panel) and Am 2 
(right panel). Unlike the case of quantum decoherence, the sensitivities to 5c achieved by the 
Kamioka-Korea setting is slightly better than those of the Korea-only and the Kamioka-only 
settings but not so much. The sensitivity is weakly correlated to 9, and the best sensitivity 
is achieved at the maximal 9. There is almost no correlation to Am 2 . 



B. Case with 5c = and 5b ^ (CPT violating) 

The allowed regions with violation of Lorentz invariance in the case of 5c = and <56 7^ 
presented in Fig. [5] have several unique features. First of all, unlike the system with 
decoherence, the sensitivity is greatest in the Kamioka-only setting, though the one by the 
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FIG. 4: The sensitivities to 5c as a function of the true (input) value of sin 2 29 = sin 2 2023 (left 
panel) and Am 2 = Am| 2 (right panel). The red solid lines are for Kamioka-Korea setting with 
each 0.27 Mton detector, while the dashed black (dotted blue) lines are for Kamioka (Korea) only 
setting with 0.54 Mton detector. The thick and the thin lines are for 99% and 90% CL, respectively. 
4 years of neutrino plus 4 years of anti-neutrino running are assumed. The other input values of 
the parameters are identical to those in Fig. O 



Kamioka-Korea setting is only slightly less by about 15 — 20%. Whereas, the sensitivity by 
the Korea-only setting is much worse, more than a factor of 2 compared to the Kamioka-only 
setting. The reason for this lies in the z/ M and z/ M survival probabilities. In this scenario, 
the effect of the non-vanishing 5b appears as the difference in the oscillation frequency 
between neutrinos and anti-neutrinos, if the energy dependence is neglected. In this case, 
the measurement at different baseline is not very important. Then the Kamioka-only setup 
turns out to be slightly better than the Kamioka-Korea setup. This case is also unique by 
having the worst sensitivity at the largest value of Am 2 (right panel). Also, the correlation 
of sensitivity to sin 2 26 (left panel) is strongest among the cases examined in this paper, 
with maximal sensitivity at maximal 6. 



C. Comparison between sensitivities of Kamioka-Korea setting and the existing 
bounds 

We summarize the results of the previous subsections in Table Ull along with the present 
bounds on 5c and 5b, respectively. We quote the current bounds on 5c's from Ref.s 38, 3^ 
which was obtained by the atmospheric neutrino data, 

|5c MT |<3xKT 26 . (11) 
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FIG. 5: The same as in Fig. 0]but for the case of non- vanishing 5b. 



We note that the current bound on 5c Mr obtained by atmospheric neutrino data is quite 
strong. The reason why the atmospheric neutrino data give much stronger limit is that the 
relevant energy is much higher (typically ~ 100 GeV) than the one we are considering (~ 1 
GeV) and the baseline is larger, as large as the Earth diameter. 
For the bound on 6b, Barger et al. [40J argue that 



|&V| < 3 x 10~ 20 GeV 



(12) 



from the analysis of the atmospheric neutrino data. Let us compare the sensitivity on 
Sb within our two detector setup with the sensitivity at a neutrino factory. Barger et 
al. [401 considered a neutrino factory with 10 19 stored muons with 20 GeV energy, and 
10 kton detector, and concluded that it can probe \5b\ < 3 x 10~ 23 GeV. The Kamioka- 
Korea two detector setup and Kamioka-only setup have five and six times better sensitivities 
compared with the neutrino factory with the assumed configuration. Of course the sensitivity 
of a neutrino factories could be improved with a larger number of stored muons and a 
larger detector. A more meaningful comparison would be possible, only when one has 
configurations for both experiments which are optimized for the purposes of each experiment. 
Still we can conclude that the two-detector setup could be powerful to probe the Lorentz 
symmetry violation. 



VI. NONSTANDARD NEUTRINO INTERACTIONS WITH MATTER 



It was suggested that neutrinos might have nonstandard neutral current interactions with 
matter 0, m,u\, 41 1, v a + / — > v@ + / (a, (3 = e, /i, r), with / being the up quarks, the down 
quarks and electrons. This effect may be described by a low energy effective Hamiltonian 
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LV parameters 


Curent bound 


Kamioka-only 


Korea-only 


Kamioka-Korea 


\8c\ 
\Sb\ (GeV) 


< 3 x 10" 26 
< 3.0 x 10" 20 


< 5 x 10" 23 

< 1 x 10" 23 


< 4 x lO" 23 
< 0.5 x 10" 23 


< 3 x 10~ 23 
< 0.6 x 10" 23 



TABLE II: Presented are the upper bounds on the velocity mixing parameter 5c and the CPT 
violating parameter 5b (in GeV) for the case where 6 m = 6^ = 9 C = 6 and rj = rj = 0. The 



current bounds are based on [38|, |40(] and are at 



/o CL. The sensitivities obtained in this study 
are also at 90% CL and correspond to the true values of the parameters Am 2 = 2.5 x 10 _3 eV 2 and 
sin 2 26*23 = 0.96. 



for new nonstandard interactions (NSI) of neutrinos: 

H mi = 2V2G F {V alp v p ) (e f j' L f L Yf' L + e f l R f R Yf R ) + H.c. 



(13) 



where = + and e^p = s{fg It is known that the presence of such NSI can affect 
production and/or detection processes of neutrinos as well as propagation of neutrinos in 
matter. In this work, for simplicity, we consider the impact of NSI only for propagation. By 
using e a p denned as e a p = J2f= u de £ a/3 n f / n e, the effects of NSI may be summarized by a 
term with dimensionless parameters e a p in the effective Hamiltonian 



-ee ^-efi ^er 

# pff = V2G F N e I e* e, 



'flfj, £[1T 



(14) 



£ £ 



£ 7 



which is to be added to the standard matter term \f2GpN e diag.(l, 0, 0) [5| in the evolution 
equation of neutrinos. Here, Gf is the Fermi constant, N e denotes the averaged electron 
number density along the neutrino trajectory in the earth. The existing constraints on e a p 



are worked out in [42|, |43|, |44( : 



-0.9 < e ee < 0.75 |e e „| < 3.8 x 10" 4 \e er \ < 



0.25 

-0.05 < £ m < 0.08 |£ WT | < 0.15 



£ 



(15) 



TT ~ 



< 0.4 



Note that the bounds on e^ T and £ 



are coming from atmospheric neutrino data [44j and 
LEP data [43f| , respectively, which are relatively weak and we wish to investigate how much 
we can improve these bounds at the Kamioka-Korea two detector setup. 

In this work we truncate the system so that we confine into the \i— r sector of the neutrino 
evolution. Then, the time evolution of the neutrinos in flavor basis can be written as 



4r" 

at \ z/. 



U 



2E 



'flT 



Ett £ 



/'/' 



(16) 



where U is the flavor mixing matrix and a = \f2GpN e . In the 2-2 element of the NSI term 
in the Hamiltonian is of the form 

&tt £ (j,fi because the oscillation probability depend upon 
e's only through this combination. The evolution equation for the anti-neutrinos are given 
by changing the signs of a and replacing U by U*. 
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In fact, one can show that the truncation to the 2x2 sub system is a good approxima- 
tion. In the full three flavor framework the disappearance oscillation probability can be 
computed to leading order of NSI as [45[ 

P(^-n> M ) = 1 - sin 2 2^ 2 3 sin 2 A 32 

r 2 
sin 2 2# 23 sin 2A 32 + cos 2 29 23 -— sin 2 A 



- \e^ T \ cos 4>fj,r sin 26> 2 3 (aL) 



- {e TT - £/i|t) sin 2 26*23 cos 26*23 (aL) 
■Am 2 



A 32 

^32 



2 2 
sin 2 A 32 - — — sin A 32 
A32 



where A 32 = A ™g L and Mr is the phase of e^ T . The result in (1171) indicates that the 
truncation is legitimate if e's are sufficiently small. Also note that the dependence on 
£tt — goes away for sin 2 ^ = 0.5, so that the muon disappearance becomes insensitive 
to this combination of e's. We will find that the sensitivity on e TT — e w strongly depends 
on sin 2 9 (9 being maximal or not) for this reason. In the following we set e w = so that 
the sensitivity contours presented for e TT actually means those for e TT — e w . Moreover, for 
simplicity, we assume that e MT is a real by ignoring its phase. 

In Fig. [61 presented are the allowed regions in e MT — e TT space for 4 years neutrino and 
4 years anti- neutrino running of the Kamioka-only (upper panels), the Korea-only (middle 
panels), and the Kamioka-Korea (bottom panels) settings. The input values e^ T and e TT are 
taken to be vanishing. 

As in the CPT-Lorentz violating case studied in Sec. IV Bl and unlike the system with 
decoherence, the Korea-only setting gives much worse sensitivity compared to the other two 
settings. Again the Kamioka-only setting has a slightly better sensitivity than the Kamioka- 
Korea setting. However we notice that the Kamioka-only setting has multiple e TT solutions 
for sin 2 #23 = 0.45. The fake solutions are nearly eliminated in the Kamioka-Korea setting. 

The sensitivities of three experimental setups at 2 a CL can be read off from Fig. [61 The 
approximate 2 a CL sensitivities of the Kamioka-Korea setup for sin 2 9 = 0.45 (sin 2 9 = 0.5) 
are: 

|e MT | < 0.03 (0.03), \e TT \ < 0.3 (1.2). (18) 

Here we neglected a barely allowed region near |e rr | =2.3, which is already excluded by the 
current data. Note that the sensitivity on e TT becomes weak for maximal mixing (sin 2 6 = 
0.5), for the above mentioned reason (see Eq. ffTTj) and the subsequent discussions). The 
Kamioka-Korea setup can improve the current bound on |e MT | by factors of ~ 5 whereas the 
bound on |e rr | we obtained is comparable to (worse than) the current bound for sin 2 9 = 
0.45 (sin 2 9 = 0.5). A similar statement applies to the case for the Kamioka-only setup. 

There are a large number of references which studied the effects of NSI and the sensitivity 
reach to NSI by the ongoing and the various future projects. We quote here only the most 
recent ones which focused on sensitivities by superbeam and reactor experiments [45| and 
neutrino factory (46[. The earlier references can be traced back through the bibliography of 
these papers. 



By combining future superbeam experiment, T2K [20j and reactor one, Double-Chooz 
23| . the authors of [45[ obtained the sensitivity of \e^ T \ to be ~ 0.25 when it is assumed to 



be real (no CP phase) while essentially no sensitivity to e TT is expected. The same authors 
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FIG. 6: The allowed regions in er„ T — e TT space for 4 years neutrino and 4 years anti-neutrino 
running. The upper, the middle, and the bottom three panels are for the Kamioka-only setting, 
the Korea-only setting, and the Kamioka-Korea setting, respectively. The left and the right panels 
are for cases with sin 2 9 = sin 2 623 = 0.45 and 0.5, respectively. The red, the yellow, and the blue 
lines indicate the allowed regions at la, 2a, and 3a CL, respectively. The input value of Am 2 2 is 
taken as 2.5 x 10~ 3 eV 2 . 



also consider the case of NOi^A experiment |21[ combined with some future upgraded reactor 



47I . |48| and obtained e„ T sensitivity 



experiment with larger detector as considered, e.g., in 
of about 0.05 which is comparable to what we obtained. 

While essentially no sensitivity of e TT is expected by superbeam, future neutrino factory 



with the so called golden channel v e 



and v e 



could reach the sensitivity to e TT at 
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the level of 
derived in 
to e 



~(IT 

obtained for e 



0.1-0.2 [46 . Despite that the sensitivity to by neutrino factory was not 
from Fig. 1 of this reference, one can naively expect that the sensitivity 
is similar to that of e ee which is ~ 0.1 or so. We conclude that the sensitivity we 
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/IT 



is not bad. 



VII. CONCLUDING REMARKS 

The Kamioka-Korea two detector system was shown to be a powerful setup for lifting 
neutrino parameter degeneracies and probing CP violation in neutrino oscillation. In this 
paper, we study sensitivities of this setup to nonstandard neutrino physics such as quantum 
decoherence, tiny violation of Lorentz symmetry, and nonstandard interactions of neutrinos 
with matter. In most cases, two detector setup is more sensitive than a single detector 
at Kamioka or Korea, except for the Lorentz violation with 5b ^ 0, and the nonstandard 
neutrino interactions with matter. The sensitivities of three experimental setups at 90% CL 
are summarized in Table [J and Table [Til for quantum decoherence and Lorentz symmetry 
violation with/without CPT symmetry, respectively. We can say that future long baseline 
experiments with two detector setup can improve the sensitivities on nonstandard neutrino 
physics in many cases. We believe that it is a useful addition to the physics capabilities of 
the Kamioka-Korea two-detector setting that are already demonstrated, namely, resolution 
of the mass hierarchy and resolving CP and the octant degeneracies. 
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